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Abstract 

We show that the Krylov-BogoUubov-Mitropolsky averaging in the canonical formulation can be used 
as a method for constructing effective Hamiltonians in the theory of strongly correlated electron 
systems. As an example, we consider the transition from the Hamiltonians of the Hubbard and 
Anderson models to the respective Hamiltonians of the t-J and Kondo models. This is a very general 
method, has several advantages over other methods, and can be used to solve a wide range of problems 
in the physics of correlated systems. 
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Model Hamiltonians used to describe strongly correlated electron systems with the electron 

potential energy much greater than its kinetic energy can be significantly simplified by reducing 

them to effective Hamiltonians in spin variables. Such a simplification can be realized by the initial 

Hamiltonian renormal-izations aimed at eliminating high-energy states and at passing to a subspace 

with lower energies of the quantum states. Eliminating high-energy states is justified in this case 

because the system properties (e.g., electrical conduction, magnetization, etc.) under the usual 

laboratory conditions are determined by the ground state and the low-energy excitations. As a rule, 

the effective Hamiltonians are constructed using the canonical transformation method [1], [2], which 

allows eliminating an "off-diagonal" small perturbation operator, which is responsible for transitions 
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between low- and high-energy states, from the original Hamiltonian. Sometimes, the perturbation 
theory is also used in systems with degenerate states [3]. Here, we propose a method for constructing 
effective Hamiltonians based on using the Krylov-Bogoliubov-Mitropolsky (KBM) averaging method 
[4] developed for solving problems in the theory of nonlinear oscillations. The averaging method 
allows eliminating the rapidly oscillating terms (the high-energy states) in any perturbation order in 
the Hamiltonian written in the interaction representation, which leads to an effective Hamiltonian 
written in the approximately "diagonal" form. In what follows, as examples, we show how the KBM 
averaging method in its canonical formulation can be used to obtain the t-J model Hamiltonian from 
the Hubbard Hamiltonian and the Kondo model Hamiltonian from the Anderson model Hamiltonian. 

The paradigmatic model in the theory of strongly correlated electron systems is the Hubbard 
model. The Hubbard Hamiltonian is written as [1], [2]: 

H = -t J2 {ctcj^+H.c.) + uY,nitnii = Ht + Hu, (1) 

where Ht is the kinetic (band) term describing the motion over the lattice sites, Hu is the operator of 
the Coulomb energy of repulsion of two electrons at the same site; is the operator of creation of an 
electron with the spin a at the site i, and t is the matrix element of electron transition from a given site 
to its neighbor. In the case of large Coulomb energy U, the appearance of two electrons at the same 
site is energetically unfavorable, and the original band splits into two Hubbard subbands: the upper 
and the lower (with a gap between them) corresponding to the one-electron and two-electron states. In 
the case of a half-filled band (one electron at the site, n = J^a- ~ Mott transition occurs for 

U ^ t, i.e., the dielectric ground state appears, and an indirect exchange coupling of antiferromagnetic 
type is established between the electrons at the site. The so-called t-J model is thus realized, which 
describes the propagation of holes in the lower subband against the background of interacting spins for 
n < 1. If the t-J model is derived from the Hubbard model, then the Coulomb term is usually taken as 
the zeroth-order approximation, and the kinetic term is considered a perturbation. In the interaction 
representation, where Hu is taken as the zeroth-order Hamiltonian, the kinetic term is divided into 
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parts describing the energetically more favorable processes (in the interior of the Hubbard subbands) 
and the less favorable processes (between the Hubbard subbands). Indeed, in this representation, 

where a = —a and the kinetic term Ht{T) = e'^^^-f/te"'^^"^ has the manifestly multiparticle 
character: 

Htir) = -t ~ ni^)4aCja{'^ - nja) + ni^cfcjallja + Hx.]- 

<ij>,<T 

-t ^ [niaCf^Cjai^ - rija) + UjaC^^Ciail - nia)Y^'' - 
<ij>,i7 

~* XI ~ nia)cl^Cjanja + (1 - nja)cj'^Ci„ni^]e~''"'^ = 
<ij>,a 

^H^ + H+{t)+H-{t). (2) 

The first and second terms in are responsible for the electron kinetics in the respective lower and 
upper Hubbard subbands; H^{t) describes the appearance of the second electron at a site already 
containing an electron, i.e., the transition from a low-energy state to a high-energy state (from the 
lower Hubbard subband into the upper); and H^{t) represents the reverse process in which the 
number of doubly occupied sites (sites with a pair of electrons) decreases by unity. The existence of 
the rapidly oscillating factors e^^^'^ in Hf^{T) and H^{t) additionally indicates that the processes 
with transitions between the Hubbard subbands are unfavorable. 

The KBM averaging method allows eliminating the rapidly oscillating "off-diagonal" terms (r) 
in any perturbation order in t (more precisely, in t/U). and can be applied to Hamiltonian (2). We 
briefly describe the method in the canonical formalism [4]- [6]. 

We consider the Liouville equation for the density matrix of a quantum system 

^ = -i[Ho + H^,p\ = -i{Lo + Li)p, (3) 

where Hq is the unperturbed "diagonal" Hamiltonian of the system. Hi is a small perturbation, i.e., 
the "off-diagonal term," and Lq and Li are the Liouvillians corresponding to the Hamiltonians. It 
follows from the condition ||i?o|| ^ ll^ill where || . . . || denotes the value of an operator in frequency 
units, that the fast motion with the period 27r/||i?i|| is superimposed on a slower process characterized 
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by time of the order of ^ ||i/i|| ^. B ([3]) We can pass to the interaction representation a — e*^°*Pi 
Li(t) e*-^«^Lie"*^o^ in Eq. (3) and then apply the KBM averaging method to the equation 

% = -^L,{r)a (4) 

in order to ehminate the rapidly oscillating terms and to construct an approximately 
"diagonal" effective Hamiltonian (Liouvillian) . We briefly describe this procedure in the simplest 
form. 

We note that because the function Li{t), is periodic, it can be expanded in the Fourier series 

Li(r) = ^Ll")e-"^ i^") = ^ f drL,{r)e-^--\ (5) 

where uj^ = 2Tin/T and T is the period, which, in particular, can coincide with the period 27r/||_ffoll or 
be multiple of it. Because we consider not the fast time "vibration" of the density matrix but its slow 
evolution, i.e., the motion averaged over several time periods T . it is natural to define the averaging 
operation 

P^a{r) = ^ 1^ dra{r)^{a{r)), (6) 

where the projection operator P'^ is the operator of averaging the rapidly varying quantities. We also 
define = 1 - (it is easy to see that P^P^ = P^, Q^Q^ = Q^, and P^Q^ = Q^P^ = 0), We can 
then write 

a{T)=P-<j{T)+Q-a{T), (7) 

i.e., we decompose the real motion of the system described by the density matrix (t(t) into the 
averaged P'^a and fast "vibration" Q'^a. Substituting expression (7) in (4) and acting on it from the 
left successively by the operators P^ and we obtain the two coupled equations 

|- (PV) = -zP^Pi (t)PV - zP^Li (t)Q V, (8) 

^{Q^a) = -iQ^Li{T)Q^a - iQ^Li(r)PV. (9) 

OT 

Equations (8) and (9) were derived using the periodicity property of the function Q'^a and the fact 
that P'^cr is a slow function of time. The solution of Eq. (9) can be represented formally as 

Q^a = -iJ dr'Q^' Li{t')P^' (t{t') ~i j dr' Q^' Li{t')Q^' <t{t'), (10) 



and the operator constant in this expression is assumed to be zero: we here use the freedom to 
choose this constant arbitrarily because one first-order differential equation (4) was divided into two 
equations, Eqs. (8) and (9), at the preceding stage. We iterate expression (10) and obtain the power 
expansion in Li: 



Q^a^!^-if dr'Q^'Liir') - i f dr' f dT"Q^' Li{t')Q^" Li{t") . . . , |(a 



), (11) 



where the slowly varying function P'^a = {a) is brought outside the integrand. Substituting formula 
(11) in (8), we obtain an equation for slowly varying quantities. For example, in the second order in 
Z/i , we have 

^(a) = -^{L,{r)){a) - Q\r' {L,{r')iL,{r') - (Li(t')))}) (a) ^ -zL,ff{a), (12) 



where we introduce the effective Liouvillian 



Leff = (Li(t)) - ^ l^f dT'{L,{T'){L,{T') - (Li(r')))}^ 



(13) 



which is associated with the effective Hamiltonian 



f dT'{m{T')-{m{r'))),m{T) 



(14) 



To derive expression (14) from (13), it is most convenient to use Fourier expansion (5). 

In the Hubbard model considered here, averaging the Hamiltonian Ht{T) over the period 2'jt/U 
results in the relation (fft(r)) = because the averages of Hf{T) are zero because of the factors 
g±ic/r ^jjg total KBM averaging procedure, i.e., formulas (14), to obtain 

= H? + ^[H+{0),H,-m. (15) 

We note that the canonical transformation method applied to the Hubbard Hamiltonian leads 
to a similar effective Hamiltonian but with additional "off-diagonal" terms ~ [H^{0), H^] (see, e.g., 
[1], [2]), which can be neglected because they take the inter band transitions into account only in the 
second order m t/U [2] but not in the first order. In the KBM method, such terms do not appear, 
because the rapidly oscillating term [H^ (r) — (r) , H°] in (15) is zero after averaging over the period 



27r/?7. Further, to obtain the desired result, it is necessary to substitute given in (2) in expression 
(15) and to perform rather simple transformations, which are described in detail in [Ij. Namely, it is 
necessary to commute the operators, to omit the three-site terms [1], [2], and to project the obtained 
Hamiltonian on the lower Hubbard subband, i.e., to omit the terms multiplied on the left and on the 
right by the operators njg- and njg such as, for example, the second term in (2): TiigC^^Cjanja (the 
high-energy motion of electrons in the upper Hubbard subband is described by similar terms). Ihe 
obtained effective Hamiltonian H^/f is just the Hamiltonian of the t- J model: 

H^ff ^ Ht-j = Eld - nia)ctcj4l - nra) + h.c.\ + JY^i^iSj - ^mrij), (16) 

where 

crcr' 

r is a vector composed of Pauli matrices, rij = rijo- and J = /U. 

The algorithm described above for constructing an effective Hamiltonian can also be used to 
transform the Hamiltonian of the Anderson model into the Hamiltonian of the Kondo model. 
Transitions from the Anderson model to the Kondo model are usually performed using the SchriefFer- 
Wolff transformation [7], [8], i.e., the canonical transformation with an anti-Hermitian operator whose 
form must be chosen. 

The KBM averaging method also has some obvious advantages in this case. The original Anderson 
Hamiltonian can be written as [7], [8] 



H^Ho + V, (17) 

Ho = Yl ^kct^Cka +sYl ^aCda + U cJ^QtcJ^Cd^, (18) 
k,tT <y 

y = Y. (^fe^'^aCda + h.c) , (19) 

where c^^ and are the creation operator and the energy of an electron with the momentum k and 
spin u in the conduction band, cj^ and e are the creation operator and the energy of a localized electron 
of the impurity atom, and U is the Coulomb interaction energy between two electrons occupying the 
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impurity atom. Here, V describes the coupling that mixes the impurity states with band states, and 
Vkd is the coupUng constant. 

To pass to the interaction representation V{t) = e^Hory^-iHoT ^ gj.g^ write Cka{T) = 
e^"°'^Ckae~^"°'' and Cdair) = e^^o'^Cd^^e"*^"'^ as 

CkAr) = Ckae-'"'^ Cdair) = Cda{l - nda)e-''' + CMae-'^'+''^^ . (20) 

As a result, we obtain 

yij) = E {^w[e'(^'=-^)"c+ Cd<,(l - m-,) + e'^'^-'-^^-c+^Cdanda] + H.c] . (21) 

k.(j 

We take all energies relative to the Fermi level. Then, in the temperature region under study, the 
energies (or £fe — instead of e^, where /z is the Fermi energy) are small, e is negative, and U, 
\e\ > Efe. 

The time average of the Hamiltonian V(t) is zero in the first-order of the perturbation theory 
because of fast oscillations of the factors e='=*(^'= and e='='(^'=~^~^)'^. In the second order, the KBM 
average of (14) is 

J\t'V{t'),V{t) ^ = 

k,k',a 

k^k' ^(T 

- ^ E ^' 0)"'^- - \ E[-^(^' ^' ^) - ^' 0)] Wda, (22) 

where 

/(fc, k', u) = VkdVC,, (- — + - — ^-—j) , 

\ek-e-U ek'-s-UJ 
and 6*^^'="^'='^'^ are slowly varying factors. Passing to the Heisenberg representation in expression (22), 
we obtain the effective Hamiltonian 

H^f^ir) ^ H^f^ ,, = Ho + ifS(O), (23) 
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where Hq is defined by (18). If the Schrieffer- Wolff transformation is used, then Hamiltonian (23) 
contains additional terms of the form [7], [8] 




which describe the higli-cncrgy processes, namely, the variation in the impurity level population 
due to the capture of two band electrons or the transition of two electrons from the d-orbital into 
the conduction band. The effective Hamiltonian in the KBM method does not contain such terms, 
because they are rapidly oscillating terms in the interaction representation and hence vanish under 
the averaging procedure. We neglect the terms describing the potential scattering of the electron 
conduction (the second term) and the electron energy renormalization on the impurity atom (the 

(2) 

third and fourth terms) in iJg^'y(O) and also the Coulomb (high-energy) term in Ho and then obtain 
the Hamiltonian of the Kondo model from (23) with expressions (18) and (22) taken into account: 

k,<T k,k' ,a 

(25) 

where the interaction term can be expressed in spin variables of the band electrons and electrons on 
the d-orbital of the impurity atom [7], [8]. The KBM averaging method can also be used to consider 
systems under nonequilibrium conditions; for example, it can be used to study the Kondo effect in the 
case where the energy of an electron localized at a quantum dot is modeled by an external alternating 
electric field. The KBM averaging is a reliable method for constructing effective Hamiltonians for 
strongly correlated electron systems: high-temperature superconductors, oxide magnetics with colossal 
magnetoresistance, quantum dots, etc. We have illustrated the use of this method with an example 
of the Hubbard Hamiltonian transformation to the Hamiltonian of the t-J model and the Anderson 
Hamiltonian transformation to the Hamiltonian of the Kondo model. The computation algorithm 
is rather simple and natural. This method does not require sophisticated tricks, which are needed 
for choosing the specific form of the unitary operator if the canonical transformation is used. It is 
also unnecessary to know the eigenfunctions and eigenvalues of the Hamiltonian in the zeroth-order 
approximation, in contrast to the case of the perturbation theory used in systems with degenerate 
states. In addition, the new method is self-sufficient: it does not result in the appearance of "off- 
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diagonal "terms of higher order in the perturbation parameter in the transformed Hamiltonian, which 
typically appear if the canonical transformation is used and are then dropped based on some plausible 
reasoning. This approach is very general and can be used to solve a wide range of problems in 
the physics of strongly correlated electron systems. In particular, it can be used to derive the 
effective Hamiltonian (of the t-J-model type) for the recently discovered class of high-temperature 
superconductors, which are iron-based layered compounds (see [9]). 
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